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q I Abstract 

Q ■ In this paper, using of the rigorous statement and rigorous proof the Maxwell 

distribution as an example, we establish estimates of the distribution depending on 
| the parameter N, the number of particles. Further, we consider the problem of 

the occurrence of dimers in a classical gas as an analog of Bose condensation and 
establish estimates of the lower level of the analog of Bose condensation. Using 
of the dequantization principles we find the relationship of this level to "capture" 
theory in the scattering problem corresponding to an interaction of the form of the 
Lennard-Jones potential. This also solves the problem of the Gibbs paradox. 

We derive the equation of state for a nonideal gas as a result of pair interactions 
£3 \ of particles in Lennard-Jones models and, for classical gases, discuss the A-transition 

to the condensed state (the state in which V^ p does not vary with increasing pressure; 
g . for heat capacity, this is the A-point). 

g : 1 The Maxwell distribution 

It is an old misconception that statistical physics and thermodynamics can be derived 
t-h ■ from the laws of mechanics and dynamical systems. It still persists from the days of the 
controversy between L. Boltzmann with H. Poincare, E. Zermelo, and other mathemati- 
cians. However, in order to solve the clusterization problem, computer simulation based, 
rule, on the laws of mechanics is usually used. For example, from a mechanical point 
of view, to obtain a dimer, i.e., a coupled pair of particles, a "collision" (interaction) of 
three particles must occur (so that part of the energy is imparted to one of them, and only 
this results in a mechanical capture). Although modern computers are powerful, one can a 
priori expect the wrong answer. The main equation of statistical physics, the Boltzmann 
equation, cannot be obtained, in principle, solely from mechanical laws. In principle, it 
^ can be derived (but not in the near future) from quantum field theory, in which there 
are no a priori prescribed pair interactions. Here we rely on the natural axiom of the 
existence of a mean field formed by N. particles, probabilistic number theory, and the 
theory of white noise at a given temperature [1] . 

Academcian N.N.Bogolyubov used to say: "I looked fora small parameter during my 
whole life." Bogolyubov was a mathematician in essence and looked for small parameters 
in physics. Physicists operate very well with digital values and intuitively understand or 
reckon up mentally whether or not a given asymptotic is applicable. 

Let us cite the corresponding text of Landau and Lifshits in their manual on statistical 
physics. Assuming that the Russel-Sounders case of connection in the atom holds, the 
authors represent the partition function in the following form (we simplify their represen- 
tation): 
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where the symbols Ej stand for the components of the fine structure of the normal term. 
Let us quote: "As is known, the existence of nuclear spin leads to the so-called hyperfine 
splitting of atomic levels. However, the intervals of this structure are so tiny that they 
can be regarded as small intervals as compared with T for all the temperatures for which 
the gas exists as gas." This is continued in a footnote: "The temperatures corresponding 
to intervals of the hyperfine structure of diverse atoms are beyond the limits from 0, 1° 
to 1,5°." ([2], Russian p. 163). Thus, the authors say that the value T = 0, IK is large, 
which just means that one must introduce the small parameter indicated above. 

This example shows that physicists do not need this parameter. Using at appropriate 
places the digits of hyperfine structure, the Russel-Sounders connection, and the digits 
arising in the specific problem under consideration, they know whether or not a given 
formula is applicable better than the mathematician who obtained the related estimates. 
However, to be correct, this is true for the Great Physicists only, and the above instruments 
can create far-reaching errors of ordinary good physicists (see, e.g., [3]). 

Let us consider a classical gas. 

The Maxwell distribution is of the form 

Ndv = [2^f) z- [m{vl+vl+vl))/2kT dv x dv y dv z , (1.2) 

where v is the velocity, T is the temperature, k is the Boltzmann constant, and Mdv is 
the relative number of particles contained in the interval dv = dv x dv y dv z . 

We must bear in mind that its interpretation as a distribution density is false in the 
general case, but valid only in the cumulative variant. This implies that the integral of the 
density (11.21) over any narrow finite velocity interval bounded below determines, indeed, 
the relative number Af VlV2 of particles in this velocity range. 

We obtain the usual Maxwell distribution for a sufficiently narrow velocity interval: 
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where N VX ^ >V2 ^ is the relative number of particles with velocities in the interval |ux|-y//c, 
|v 2 |\/k, and k is a small parameter. 

By a change of variables, we can get rid of the small parameter k by transferring it to 
the exponential, i.e., 

f lv2{ e-( Kmv2 V 2kT v 2 d\v\ 

^1*2 = J~ (Km v*)/2kT v 2 d \ v \ > fOT K ^ - ( L4 ) 

This formula, as an asymptotic formula, will be obtained below, as well as its estimate, 
i.e., its domain of applicability. 

Consider the most often used Lennard- Jones interaction potential 



<p{r) = 4C/ 
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where a is the distance at which the potential function changes sign and Uq is the minimum 
value of the potential (at the point r = 2 1 / 6 (j) or the depth of the potential well. 

From dimensional considerations fll.2j) for the quantities appearing in the definition of 
the particle (neutral molecule), we can write 

\ 3/2 /»oo 

^ / e-^ m ^' 2kT dv x dv y dv z = N. (1.6) 



for this molecule. Hence 

4 ir (2kT\ 3/2 
3 ' ^2 \U^) 



N, (1.7) 



Therefore, the velocity interval on which we can determine the relative number of particles 
is 

where V\ and t> 2 arbitrary velocities, |t> 2 | > |t>i|, independent of the number N. 

For a rigorous justification of the Maxwell distribution, we use the seemingly insignif- 
icant fact that the number of particles iV is an integer and apply number theory, which 
does not seem relevant at all. 

The Maxwell distribution is equally important in complexity theory as the Poisson, 
Gauss, and other classical distributions. 

Let us define the energy 

47T/3 ■ (kT/2f/ 2 

E = ~m • (L9) 

^0 

The distribution in the energy interval between mvf/2 and mv 2 /2 must have the following 
form: 



(1.10) 



Without loss of generality, let 

^ = Wo, = + !)«,, 

where / is an integer. Hence, up to 0(k 2 ), we can write 

= u^m + 0{k2)iV2 Uo (1 n) 

EE E 

If we split the total energy E into intervals of the form {IUq, (I + l)?7o}> I — 0,1, ... ,Ie, 
so that the sum of these intervals is less than E by a quantity O(k), then, in view of the 
Euler-Maclaurin formula, we have the order 

l E = K~ 7 I\ 

Further, if we replace I 1 / 2 by its integer part [I 1 / 2 ], then we decrease the sum of the 
intervals by at most a quantity 0(/t~ 7//5 ). 

Hence the union of the partitions Uq[1 1 I 2 }1 satisfies the inequalities 

Ie 

E-O(k' 7 / 5 )<kU J2[ i1/2 ) 1 < e - (1-12) 

o 

Thus, we obtain energy boxes and wish to find the most probable number of particles 
with energies in each box. 



Now let us split the number of particles N = ^Njk, where k — 1, 2, . . . , [j 1 ' 2 ), j = 
1;2, ...,Ie- 

Therefore, given condition fl!.12j> . we obtain the following constraint on our partition: 

l E li 1/2 ] 

E - 0(kT 7/5 ) < kU N i k - E - ( L13 ) 

j=l k=l 

The condition N = ^jk implies that the size of the ordered sample with replace- 
ment |4] is equal to N, while condition f II . 1 3 j) means that the energy corresponding to 
this sample is contained in the interval {E — 0(k~ 7 ^ 5 ), E}. 

These were heuristic considerations. Now we make the following assumptions. 

In the volume V, consider the system of N particles possessing the energy E. Moreover, 
N -> oo, E/Uq -> oo. 

The interval (0,E) is divided into small (compared with E) subintervals E^ < E i+ i, 
i = 1, . . . , lo, and the corresponding intervals of the moduli of velocities \vi\ < fj+i, as well 
as intervals of the phase volume, and the energy boxes AQi that are contained between 
these velocities 

AQj = f dv x dv 2 dv z = const(Eff 1 - E- /2 )V sa constA£<£j /2 V. (1.14) 

J\Vi\ 2 

In these energy boxes, we place different particles using all possible ways N. In other 
words, we take an ordered sample with replacement from iV "balls" to these energy boxes 
(phase volumes) by the method indicated in (11.131) : 

lE [i 1/2 l 

< kU Yl N i k - E - ( L15 ) 

3=1 k=l 

By A/"(|fi|, [^2 1 ) we denote the relative number of particles in the velocity interval 

(ui, i^- 
Under the conditions given above, the following theorem is valid. 

Theorem 1.1. The probability that the estimate 

rM + 0(N-V^) 2/2kTd 

11 ' 1 ^ J-°°e-™ 2 / 2kT dv 1 dv 2 dv 3 



0[ ^-^|ln In N\ £ ), (1.16) 



N 

where 5 > 0, and e is any arbitrarily small number, and \vi\ > are arbitrary velocities), 
does not hold is exponentially small (is less than 1/N k , where k is any integer). 

The same assertion is also valid for any large velocity interval, i.e., v 2 > 5 > 0, where 
S is independent of N. 

By analogy with the term "convergence in measure," we can state that, in (11.161) . 
there is an "estimate in measure." In fact, the physical formula (ll.2p can be rewritten in 
the more exact form 

r v 1 +0(N- 1 / 2 + s ) / \ 3/2 

N Av = J \^r) ^ {m{vl+vl+vl))/2kT dv x dv v dv z , (1.17) 

where A = 0(A^ _1 / 2+5 ) and 5 > 0. 

In usual probability notation, the theorem can be restated as follows. 
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Theorem 1.2. The following relation holds: 

/ rW\+o{N^)+ s ) 2 /2kT 



E 



f-™e-™ 2 / 2kT dv l dv 2 dv 3 ~ VN 

= 0(N- k ), (1.18) 

where k is any number, e > is an arbitrarily small number, 5 is any number, and \v±\ > 
are arbitrary velocities. Here P is the Lebesgue measure of the phase volume defined in 
parentheses in (jl.lBp with respect to the total volume. 

These estimates are sharp (unimprovable). The theorem belongs to number theory. 
It has no relation to particle dynamics in which the Maxwell distribution is derived from 
the Boltzmann equation, which has not been is rigorously justified up to now. The usual 
dynamical approach and its criticism is contained in Kozlov's book [3]. 

Nevertheless, it is natural that, under certain conditions, the dynamical system attains 
the most probable (from the point of view of probability number theory) distribution. This 
consideration can be useful for the dynamical approach. 

Let us present a sufficiently elementary proof the theorem on the Maxwell distribution, 
without, essentially, referring to important and elegant results of number theory based on 
the Meinardus theorem [6], Theorem 6.2 and on Vershik's elegant theory of multiplicative 
measures [7], which could help us avoid some inessential and deliberate manifestations of 
integrality (such as taking the integer part of Z 1 / 2 ). 

On the other hand, the given estimates, which the author used in his papers dealing 
with economics and linguistics [8], [9], are more understandable to readers that are not 
experts in number theory and probability theory, in particular, to physicists and specialists 
in analysis. 

The proof is based on the estimates given by the author in [10] and on a theorem 
similar to the Meinardus theorem. 

For the 6-dimensional detailed proof was given in [TI] . Essentially, we repeat 

this proof for the 3-dimensional case. 

Let us study the system defined as follows. For energy levels j — 0, 1, 2, . . . of multi- 
plicities 

q l =\3 1/ % j = 0,1,2,..., (1.19) 

we consider all possible collections {Nj k } of nonnegative integers Nj k , j = 0,1,2,..., 
k = 1, . . . , qi satisfying the conditions 



Y,N 3k = N; (1.20) 



3=0 

oo qi oo qi 

j=0 k=l 3=1 k=l 

where N and M are given positive numbers (which can be assumed integers without loss 
of generality) All such collections are assumed equiprobable. 

Denote the dimensionless quantity M = E/Uq and define the numbers (3 and £ as 
solutions of the system of equations 

oo oo 

m = r 1 E-^V^ n = r 1 E^V*'. 

j=l j=0 



Remark 1. By the Euler-Maclaurin formula, we have 



poo I 

N = r 1 / x 1/2 e-^ dx{l + 0{f3)) = =- v^F/r 3/2 C\ 
Jo 2 

/■OO o 

m = r 1 / ^ 3/ v^ rfx (i + o(/?)) = - v^F/^-^r 1 - 

Jo 4 

Hence (3 and £, 

2iV _ 1 M 3 / 2 

Remark 2. As an example, we consider the Maxwell distribution. Note that the same 
argument can be used for any arbitrary Gibbs distribution, but with more cumbersome 
estimates. Just the same estimates are obtained for the Gibbs distribution corresponding 
to the Hamiltonian E = (p 2 + q 2 ) 2 , because the phase cells for this Hamiltonian satisfy 
the same relations fll.7p - fll.15p . 

We consider sufficiently general classical Hamiltonian function H(p, q), where g G i? 3 , 
p G -R 3 , under the following two assumptions: 

1) H(p, q) — > oo as \p\ + \q\ — > oo not slower than (|p| + \q\) a for some a > 0; 

2) the function 

7(A) = / dpdq= [ dE 1 1 6{E-H(p,q))dpdq, 

JH{p,q)<k JO J J 

under the assumption that V (A) is a sufficiently smooth function, determines a phase 
cell invariant under the Hamiltonian system corresponding to the Hamiltonian H(p,q). 
We choose a partition such that Ei + \ — E\ = Eq. Then 

£, = £b(Z + l). (1-22) 

Let Ni be an ordered sample with replacement to the cell £7+1 — E\. An ordered sample 
with replacement from N balls to cells invariant under the Hamiltonian system (to "energy 
boxes" ) 

"E l+1 r E o( l + 1 ) 



/ H(p,q)dpdq= / \V\X)d\ 
J Ei JE l 



leads to the state 

J2NiE lQl <£ N , gi = C[V'(E l)], (1.23) 
1 

where C is a constant. Then the proof and the estimates are just the same as in Theo- 
rem [L2J 

Denote by N{M, N) the total number of collections {A^} satisfying the con- 
straints fJOQD (OTj) . 

We assume everywhere that the parameters (3 and £ satisfy the relation 

^ < £-3/2+ £) (1.24) 

for an arbitrary (but fixed) e > 0. 
Let 

fj, = ln£, so that £ = e M . 
Suppose that A1 is the set of ordered samples satisfying conditions (11.201) and (11.211) . 



For the numbers J\f(N, M) of such variants, we obtain the following estimate: 

M(M,N) < C\fNexp{N\nN + (3M + fiN}. (1.25) 

Indeed, 

Af(MN) = W T \ < W Y eXP{ ~ Ej "° Etl NjkW + ^ 

oo <?j oo _7V- fe (/3j+/x) 00 9 -? 

= N! yiyi j2 "—^ — = m e/3M+ ^ n n 

j=0 k=l N jk =0 i k ' j=0 k=l 

= N\ exp j/3M + fiN + Y^ Qj e ~ Pj ~^\ = N - exp{/3M + fiN + N} 
< j=0 > 

< CVN exp{iVlniV + (3M + fiN} (by Stirling's formula). 

Suppose that .Ma C M. is the subset of variants such that 

i % 



EE(%-^ 

j=0 k=l 



> A, (1.26) 



where 

N ik = e-V-* = ^. (1.27) 

For the number J\f(M, N, A) of the sample from M. a, we obtain the estimate 
Af(M, N, A) < N\ exp j/3M + /iiV - cA + ^ qje~ fii ~" j 

exp|^(g J e- /3 ^^ +c - cnA + exp j^( 9i e-^- c + cJV^) 

v ^ j=Q ' ^ j=0 



:i.28i 



for < c < /i, where iVj is given by (jl.27p . 

Further, as in [121 . we take two terms of the expansion in the Taylor series 



q je - pj -^ ±c =F cNj = q 3 e- p3 -^{e ±c T c) = {\ + ^e ±0c ^j , (1.29) 



where 9{i) is some midpoint, 9 = 9(c) G (0, 1). 

If c < min{yu/2, 1}, then this implies the inequality 



i 

q^^- 6 ^ < 2Ke-"(3- 3/2 , (1.30) 

3=0 



where K is a constant. 
Therefore, 



W(iV,M,A) < CVN exp{N\nN + PM + fxN} exp{~cA + Kc^e-^p-^ 2 }. (1.31) 
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We substitute 

A = VN lniV|ln lniV| e x e'^p-^VhJj \\n\n N\ £ (1.32) 

and 



c = - - x /? 3 /V /2 VmlV|mln/V| £ (1.33) 
2K 

in (11.3ip . This implies that, in particular, 

W(iV,M,A) < C k VN exp{N\nN + (3M + fiN}N- k (1.34) 

for any k. 

Let us now find a lower bound for these quantities. 

We estimate the number of samples A/q(M, N) < J\f(M, N) satisfying conditions (ll.20p 
and (11.211) ; moreover, in the last inequality, we consider the equality 

OO 1j OO 1j 

E = EE^ = M - (i.35) 

j=0 k=l j=0 k=l 

Suppose that M.q is the set of collections of occupation numbers satisfying (11.351) . Then 

TV! 



AT (M,N)= Yl 



{N jh }eM nj=o rifc=i N jk \ 

m \- m Er=o Eti N jk )s(M, et=q gtij^) n _ 

~ 2^, rr°° IT 3 iV-J ' 1 J 

Here the sum in the second row is taken over all finite collections of nonnegative occupation 
numbers and 5(m,n) = 5 mn is the Kronecker delta. 
Substitute the integral representation 

6 mn = ^- f ' e {lsx+uj){m - n) dx 



J mn ^ 
111 



■n/s 



of the Kronecker symbol (where s and to are arbitrary nonzero real numbers) into (|L~c 
choosing s — 1 and u = fi for the first factor and s = uj = (3 for the second factor. Then, 
for A/o(M, N), we obtain the integral representation 

M>(M, N) = / / e A *^ # (1-37) 



where 

A = /r^V x TV, (1.38) 

OO 

$(y>, V>) = i(3 5/2 e>*M<p + i(3 3/2 e^N^j + (3 3/2 ^ qje -M-W+Pw) . (1.39) 

3=0 

(The sign x TV means that there exist constants c\ and C2 such that c±N < A < C2N). 



s 



Indeed, the substitution described above yields 



M~o(M,N) 



(3N\ 


e /3M+tiN 




A 9 


f3N\ 


e /3M+nN 




47T 2 


(5N\ 


e f3M+iiN 




47T 2 


f3N\ 


e /3M+nN 



E 



Jf3ipM+iipN 



j V| , i} 11/ II II/' I V /'' ! 

oo <?j 



X 



{CO ^ x 

7=0 fc=l ' ' 



Jf3tpM+iipN 



—ft/ft \J —it 



J 

oo Qj oo 

rins 

i=0 k=lN jk =0 
co 9j 



e ^ 



N jk (f3j+fj,+i(3jip+i^) 



N jk \ 



dip J cZy 



j=0 fc=l 



4vr 2 



exp<^ ifiipM + iipN ^ qj e~ w+ ^ + 
\> j=0 



if3jip+iijj) 



dip ) dip. 
1.40) 



Lemma 1. The phase function Q((p,ip) defined by (II. 39!) possesses the following proper- 
ties: 

1. All of its derivatives are uniformly bounded for the values of (3 and £ satisfying 
inequality (11.241) . 

2. The phase function has a stationary point ip = mod27r/6, ip = mod27r. 

3. The matrix $"(0,0) of second derivatives of the phase function at the stationary 
point is nondegenerate and is strictly negative definite uniformly in the parameters (3 
and £ satisfying inequality fll .24j) . 

4. the imaginary part of the phase function at the stationary point is zero and its real 
part attains an absolute maximum there; moreover, for any 7 > ; there exists 
a 5 > independent of the parameters j3 and £ satisfying inequality (ll.24p such that 



Re$((p,ip) < Re $(0,0) -5 for dist((y>, ip), (0, 0)) > 7. 



'i.4r 



Proof. 1. The boundedness of the derivatives of the phase function is proved by direct 
calculations. 

2. To to verify that the point (0, 0) is a stationary point of the phase function, let us 
calculate its first derivatives: 



d<5> 

dip 



M - e-^^jq^-^-^^ 

/3j-i(i/j+(3j(p) 



3=0 
oo 



N-e-^qe 

j=0 



1.42) 
1.43) 



For ip = ip — 0, both derivatives vanish by the definition of the parameters (3 and £ = e M . 
3. The matrix $"(0, 0) is of the form 



$"(o,o) = -^/? 3 /V-^ 

3=0 



1 Pj 

Pj P 2 j 2 



1.44) 



Let us estimate this matrix as the matrix of the corresponding quadratic form as follows: 

(1.45) 



$"(0,0) < - P 3/2 Q3 e ~ Pj 

j=[xi/P] 



where X2 > X\ > are arbitrary fixed numbers. 

For small f3, in view of the asymptotics qj ~ j/2 for large j, the matrix on the right- 
hand side can be calculated by the Euler-Maclaurin formula, obtaining as a result, up 
to o(l), the matrix 



\J Xl 



xe x dx 

x 3/2 e -x j x 



x 3/2 e~ x dx^ 



j-i 



x 2 e x dx 



i rax 



[l,x) 
2 \(x, 1) (x, x) 



1.46) 



where 



[u,v] 



■1-2 



u(x)v(x)e x xdx 



j-i 



is the inner product L 2 ([xi, X2], e~ x x). 

Since the functions 1 and x are linearly independent, the matrix (11.461) is negative 
definite, which proves the required assertion. 

4. For $>(ip,if)), from formula (11.391) we obtain 

00 

Re $(0, 0) - Re <b(tp, if>) = (3 3/2 ^ ^"^'(l - cos(^ + /3j<p)). 

j=0 

All the summands on the right-hand side are nonnegative. Therefore, omitting part of 
them and estimating the coefficients /3qje~^ for the remaining summands, we obtain 

[«a/)9] 

Re $(0,0) - Re$(<p,ip) > const /5 ^ (1 - cos(j(p + ip)) 

3 = [*l/P] 



> const — X\ — (3 — f3 



sin — 
2 



-1 



1.47) 



Choosing x\ and X2 in a suitable way, we obtain the required assertion. The lemma is 
proved. □ 

Using this lemma, we can calculate the integral ( 11.371) by the saddle-point method and 
obtain a lower bound for the number of ordered samples in the form 



N(N, M) > Cphr x VN exp{iV "In N + (5M + fxN}. 



1.48) 



Now we can estimate the integral (11.40p . By Lemma CD, all the derivative of the 
function $(<£>) are uniformly bounded. In addition, if $ is expressed in the form $ = 
$1 + i$2, where $1 and $2 are real, then, 



$;(o) = o, $1(0) < -c < 0, $' 2 (o) = $2(0) = 0. 



in 



Hence, for \<p\ < e, where e > is sufficiently small, using the Taylor formula with 
remainder, we obtain the estimates 

$(0) - CiM 2 < $i(¥>) < $(0) - C 2 \if\ 2 , (1.49) 
|$ 2 (</?)l <C 3 \<p\ 3 (1.50) 

where the Cj are positive constants independent of M and the sequence {lRjl}. 
Suppose that 

1 = 1px((p) +1p2((f) 

is a nonnegative smooth partition of unity on the circle S 1 3 ip of radius b such that 



supp^i C [—£,£] and ipi(<p) = 1 for tp G 
Let us express the integral 

1=1 exp{/r 3/2 $(^)} dip 
Js 1 



e e 
2' 2 



as the sum 



/= / exp{^ 3 / 2 <l>(^)}^iM^+ / exp{(3~ 3/2 ^)}Mv)d^ = h + h. 
Js 1 Js 1 

By Lemma [U item 3, 

Re$(y?) < $(0) -5, 5>0, 

on the support of the integrand in I 2 , while the measure of the support is of the order 
of Therefore, 

\h\ <^exp|^^ + /5- 3/2 $(0)}, b^O, (1.51) 

where is a constant. 

Let us now estimate the integral I\. For convenience, denote provisionally by h — (3 3 ' 2 
the small parameter in the exponential of our integral. On the interval D = [— e, e], we 
distinguish two subintervals D 1 / 2 C D 1 / 3 C D by setting 

D 1/2 = [-eh 1 / 2 , eh 1 / 2 ], D 1/3 = [-eh 1 ' 3 , eh 1 ' 3 }. (1.52) 

Then 

< C 3 e 3 for ip G Di /3 , 



so that (for a sufficiently small e) the imaginary part of the argument of the exponen- 
tial D\/3 is small and the following relation holds: 

Ree^' h > Ie* l(v)/ \ <p G D 1/3 . (1.53) 



Further, 



m>*M>m.c^, V e Dl/2 . (i.54) 

h h h 



1 1 



Combining this with the previous inequality and taking into account the fact that the 
length of the interval D1/2 is equal to leb}^ \ we obtain 

Re / e*M' h M<P) d V > C^' h h l l 2 . 
Jd 1/2 

Further, 

Re / e^vWfafa) d<p > 



by virtue of f 11 . 5 3 j) . Moreover, the following inequality holds: 
so that 

' < c 5 e^ h e- c ^ h ~ 1/3 



e Hv)/h ip 1 (<p) dip 



'D\D 1/3 

Combining all the previous estimates, we obtain 

I = Re I > C 6 /3 3/4 exp{/r 3/2 $(0)}. 

It remains to substitute this estimate into formula (ll.40p for J\f (M, N) and, in view of 
the formulas (jl.39j) for the phase function and by the inequality Af(M, N) > J\f (M, N), 
we obtain a lower bound for J\f{M,N). As a result, we obtain Theorem 11.11 Since the 
number M{M, N) corresponds to the Lebesgue measure of the total phase volume and the 
number J\f{M, N, A) corresponds to the Lebesgue measure of the phase volume defined 
in parentheses in formula (11.181) . we obtain the proof of Theorem 11.21 

Remark 3. The Maxwell distribution (11.21) . (11.31) holds for a "classical ideal gas" in com- 
mon understanding. By definition of pressure P of specific volume V sp = V/N for a 
"classical ideal gas" the compressibility factor 

is identically equal to 1. 



2 Clusterization in an ideal gas and dependence of 
the compressibility factor on the pressure 

Each scientist who refutes a century old theory runs the risk of being accused of incompe- 
tence and of irritating those scientists who absorbed the old theory "with their mother's 
milk." And if this is a scientist who has achieved a good deal in his area of knowledge, he 
also runs the risk of losing his hard-earned authority. This is borne out by the history of 
new discoveries in physics. Thus, the great physicist Boltzmann, virulently attacked by 
his contemporaries, committed suicide by throwing himself down the well of a staircase. 

In 1900, Planck proposed his famous formula describing black body radiation, which 
gave results coinciding with experiments, but which he had not rigorously established. The 
mathematician Bose from India noticed that, in order to derive the formula, one must 
use a new statistic instead of the old one, the so-called Boltzmann or Gibbs statistic. It 
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is possible that Planck was also aware of this statistic, but was afraid of being criticized 
or did not really believe in his own result. Bose, just like Boltzmann, was the object 
of virulent criticism, until Einstein gave his approval to the proposed statistic, which 
was also justified by the philosophical concepts of Ernst Mach. At first, physicists were 
bewildered and could not understand the Bose statistic, because they could not imagine 
how moving particles can exchange positions without using up any energy. 

These two statistics have been illustrated above by a simple financial example. 

The reply to the bewilderment of physicists was given by Mach's philosophical con- 
ception, claiming that the basic notions of classical physics (space, time, motion) are 
subjective in origin, and the external world is merely the sum of our feelings, and the goal 
of science is to describe these feelings. Therefore, if we are unable to distinguish particles 
in our subjective perception, then they are undistinguishable. 

I propose a completely different philosophy. We can regard particles as distinguishable 
as well as undistinguishable. This only depends on the aspect of the system of particles 
that we are interested in, i.e., depends on the question we are seeking an answer to. Thus, 
returning to the money example, people are interested in the denominations of the bank 
notes they own, not in their serial numbers (unless, of course, they believe in "lucky 
numbers"). 

The situation in physics is similar. Suppose we have a receptacle filled with gas 
consisting of numerous moving particles. If we take a slow snapshot of the gas, the 
moving particles will display "tails" whose lengths depend on the velocity of the particle: 
the faster the motion, the longer the tail. Using such a photograph, we can determine 
the number of particles that move within a given interval of velocities. And we don't 
care where which individual particle is located and which particular particle has the given 
velocity. 

I have derived formulas which show how the number of particles is distributed with 
respect to velocity, for example, they show for what number (numerical interval) it is 
most probable to meet a particle moving with a velocity in that interval. 

These formulas lead to a surprising mathematical fact: there exists a certain maximal 
number of particles after which the formulas must be drastically modified. If the number 
of particles is much less than this maximal number, the formulas coincide with the Gibbs 
distribution up to multiplication by a constant. Nevertheless, this is essential, because 
the corrected Gibbs formula thus obtained no longer leads to the Gibbs paradox. 

The paradox now bearing his name was stated by Gibbs in his paper " On the equilib- 
rium of heterogenous matter," published in several installments in 1876-1879, and resulted 
in great interest on the part of physicists, mathematicians, and philosophers. This prob- 
lem was studied by H. Poincare, G. Lorentz, J. Van-der-Waals, V. Nernst, M. Planck, E. 
Fermi, A. Einstein, J. von Neumann, E. Schrodinger, I. E. Tamm, P. V. Bridgeman, L. 
Brillouin, A. Lande and others, among them nine Nobel Prize laureates. 

From my point of view, the solution of the Gibbs paradox can be obtained once we 
realize that the Gibbs formula in its classical form is invalid and we modify it in the way 
that I have indicated. This modification was previously interpreted as a consequence of 
quantum theory, but this is erroneous from the mathematical point of view, since the 
passage from quantum mechanics to classical mechanics cannot change symmetry and 
therefore cannot change the statistics. 

In this situation, the following phenomenon, rather strange from the mathematical 
point of view, arises. If the number of particles is greater than the maximal number 
indicated above, then the "superfluous" particles, as we already explained, do not fit into 
the obtained distribution and the velocity of these particles turns out to be much less 
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than the mean velocity of particles in the gas. This effect differs from the Bose-Einstein 
condensate phenomenon from quantum theory, because in quantum theory these particles 
are at the very lowest energy level, they have the lowest speed, i.e., roughly speaking, they 
stop. 

Further, I try to give a physical interpretation to the obtained rigorous mathematical 
formulas. I interpret the maximal number of particles mentioned above as oversaturated 
vapor; the superfluous particles are then regarded as nuclei around which droplets begin 
to grow. As a result, this can explain the so-called phase transition of the first kind, in 
which, as the result of the system achieving equilibrium, the number of particles changes 
from that number for an oversaturated gas to that for a saturated one. Indeed, it is only 
those particles which move at speeds greater than the speed of the " superfluous" particles 
that can be doubtlessly regarded as particles of the "pure" gas (vapor), while the others 
have condensed or have mixed with the condensed particles (clusters). 

In section 1, 1 cited an example from economics, similar to the one above, that supports 
exchangeability theory (instead of the "independence condition"). In my opinion, we must 
revise, in this vein, the "Gibbs conjecture on thermodynamic equilibrium," which is based 
(see [5]) on the property of independence leading to the theorem on the multiplication of 
probabilities. It is this conjecture that leads to the Gibbs distribution, which is refuted by 
the Gibbs paradox, i.e., in essence, by the mathematical counterexample to this conjecture, 
as mentioned above. 

It is difficult for physicists to grasp this problem, because it involves a mathematical 
effect of the type of Bose condensation, which results in the appearance of a " Bose con- 
densate," which, from the author's point of view, has been treated as some coagulation 
of particles with low velocities and the formation of dimers, trimers, and other clusters. 

The phenomenon of the appearance of dimers is usually obtained by modeling involving 
the initial conditions and interactions, for example, of Lennard- Jones type. According to 
the author's point of view, if this phenomenon involves interaction, then it can occur 
before the switching-on of an interaction of Lennard- Jones type: as far as the specific 
volume is concerned, we still deal with an ideal gas. Such type of interaction is observed, 
for example, in the gas Cgo (fullerene) possessing very weak attraction (of order 0(l/r 9 )). 
It is related to the asymmetry of the molecules and the types of adjoining faces of the 
molecules of C 60 . 

This is much easier to observe experimentally, because fullerene has no liquid phase 
and is immediately transformed into fullerite particles. 

The presence of such a "saturated" total number of particles in the problem under 
consideration, with surplus particles going somewhere (passing into the Bose condensate 
0), is a mathematical fact rigorously proved together with clear estimates of where such 
aggregates may occur. However, it is not quite correct to say that the particles are 
added. Indeed, it is better to say that we lowered the temperature, while using a piston 
to maintain a constant pressure, and hence the saturated total number of particles is 
decreased. And we can simply say that, for a given temperature, the pressure is increased 
until the A-transition occurs in the "Bose condensate." The question is: Where have 
the other particles gone if the temperature is lowered simultaneously with the pressure 
limitation or the pressure at the given temperature becomes sufficiently large? Perhaps, 
they precipitate on the walls of the vessel? Such a law of "necessary" precipitation 
(coagulation) on the walls would be more interesting still and would have important 

1 The physicists to whom I described this theory warned me not to use the term "Bose condensate," 
because this evokes associations obscuring the understanding of the proposed theory. 
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practical applications. However, experiments tend to support, to a greater extent, the 
first point of view. The physicists are even of the opinion that the transition to dimers is 
a phase transition. 

The most significant fact is that this estimate not improvable. This fact follows from 
Theorem 2 in [13] . In the case of saturation, it makes it possible to determine the number 
of particles passing into clusters as the temperature is lowered, while a constant pressure 
is maintained by a piston (see [15J). This also solves the Gibbs paradox. 

The use of an unordered sample with replacement leads us to a mathematical formula 
for the Bose gas, however, without the parameter h, the Planck constant, but with the 
same parameters that appeared when using the parameters of the Lennard- Jones interac- 
tion potential. Instead of formula (11.171) . we thus obtain 



-Vi+Av / \ 3/2 

A/a^ = / I ] A ttv 2 dv. (2.1) 
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However, Av now depends on V\ in the following way: if v\ ~ 1, then Av = N Q l ^ 2+& , 
where 5 > 0, and No is the number of particles saturating the volume V at temperature T 
and maximal energy E/Uq; namely, 

at _ f m \ 3/2 r 

N °-{~u ) J ^WH' (2 - 2) 

In view of the given parameters, the velocity can be expressed as v — \JU$/2m. 
Suppose that vq is the minimal velocity; it is equal to vq = vN" 1 /'^ 5 , where 5 > 
determines the smallness of vq. For 5 = 1/3, we obtain vq = v; therefore, we set 1/3 > 
5 > 0. The estimate of the error in the formula for the distribution (12.11) is of the form 

O ( AT-V3-5/2 v /j^y Q n ln N y^ . 

namely, the following theorem is valid. 
Theorem 2.1. The following relation holds: 

I r\v \+O(N- 1/3+S i) / m \3/2 1 

P I A/ah / I — 77^ I 5-7777^ dv x dv v dv z 
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> N~ 1/3 - 5/2 Vh[N(\n\nNy 

< O(N ~ k ), (2.3) 

where k is any integer, e > is arbitrarily small, Vq > 0, < 5 < 1/3, 8% > 5, v = 
^/U /2m, and Av = 0(N~ 1/3+Sl . Here P is is the Lebesgue measure of the phase volume 
defined in parentheses in (12.31) with respect to the total volume. 

The proof of Theorem 12.11 is similar to the proof of Theorem 11.21 except that we use 
unordered samples of "balls" with replacement. 

Suppose that there is a sequence of boxes Uj, j = 0,1,2,..., and each box Uj is 
divided into qj compartments. We take N identical balls and put them into the boxes at 
random observing the only condition that 

00 

J2jNj<M, (2.4) 
3=0 

IF) 



where Nj is the number of balls in the box Uj and M is a positive integer specified in 
advance. As an outcome, we obtain a sequence of nonnegative integers Nj, j = 0, 1, 2, . . . , 
such that 



oo 

Z>i = # (2-5) 

3=0 



and condition ( 12.41) is satisfied. It is easily seen that, given M and N, there are finitely 
many such sequences. Suppose that all allocations of balls to compartments are equiprob- 
able. Since the number of ways to distribute Nj indistinguishable balls over qj compart- 
ments is equal to 

'q j + Nj-l\ T(q 3 + Nj) 

Nj J F(Nj + l)T(qj) { ■ > 

(where T(x) is the Euler gamma function), it follows that each sequence {Nj} can be 
realized in f({Nj}) ways, where 

and the probability of this sequence is equal to f({Nj}) divided by the sum of the ex- 
pressions similar to (12.7ft over all sequences of nonnegative integers satisfying the con- 
straints (12.41) and (12.51) . This makes the set of all such sequences a probability space; 
the corresponding probabilities will be denoted by P(-). The numbers qj are called the 
multiplicities. We shall assume that go is some positive integer and 

9i = b' V2 ] 5 J = 1,2,..., (2.8) 

where the brackets stand for the integer part of a number. 

What happens as M, N — > oo? It turns out that the so-called condensation phe- 
nomenon occurs: if N tends to infinity too rapidly, namely, if N exceeds some threshold 
N cr = N cr (M), then a majority of the excessive N — N cr balls end up landing in the box 
Uq\ more precisely, with probability asymptotically equal to 1, the number of balls in Uq 
is close to iV — N cr (and accordingly, the total number of balls in all the other boxes is 
close to N cr , now matter how large N itself is). Let us give the scheme of proof analogous 
to the proof of Theorem 11.21 

Define N cr = N cr (M) by the formula 
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where b is the unique positive root of the equation 



J2J^ = M. (2.10) 

Next, let 
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A = iV c 2 / 3+£ , (2.11) 

where e > is arbitrarily small (but fixed). If N > N cr , then there exist constants C m 
such that 

P{\N -{N-N cr )\ > A) <C m N~ m , m = l,2,... . (2.12) 
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It is not hard to compute N cr (M). Indeed, in view of (I2.8p . the Euler-Maclaurin 
formula gives 

Y.Jrbi- f ^ = r^a\) (2.13) 

• 1 J 



(where ((x) is the Euler zeta function) and likewise 
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Eprrr-fl- 3 ' 2 / ^rff = 3- 3/2 r(i)C(i). (2.14) 

By substituting this into (12.91) and (12.101) . we obtain 

A^r(|) C (l) 

N " mmw ( 5) 

In contrast to the Maxwell distribution, the compressibility factor for the given distri- 
bution equals 

f (p 2 /2m)p 2 dp 

Z=^r = ~ J ^Z- 1 = 0.523. (2.16) 
However, if the number of particles Af C ^0, then 

(p 2 /2m)p 2 dp 
e (p2/ 2m _ M )/fcT_ 1 



py 2 ./' — 
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(2.17) 



(p2 /2m-/j)/fcT _J 



i.e., there appears a negative parameter /z which tends to —00 (and Z — > 1) as N/N Q 
decreases to zero. 

Consider the gas which consists of different molecules, or dimers, trimers, . . . , 
fc-mers. 

Now suppose that the situation is the same, but we should additionally paint each of 
the N balls at random into one of K distinct colors. Now that we can distinguish between 
balls of different colors but balls of a same color are indistinguishable, how does this affect 
the probabilities? 

Instead of immediately painting the balls, we can further divide each of the qj com- 
partments in the jth box into K sub-compartments and put the uncolored balls there 
(with the understanding that the balls in the kth sub-compartment will then be painted 
into the kth color and the dividing walls between the sub-compartments will be removed). 
Now we have Kqj sub-compartments in the jth box, so that there are 

Kq j + N j -1\_ TiKq. + N,) 



N 5 J V(N j + i)r(K qj ) 

ways to put Nj balls into the jth box. All in all, the introduction of K colors has the 
only effect that all multiplicities qj are multiplied by K. 

Our theorem applies in the new situation (with qj replaced by the new multiplicities 
qj = Kqj). The computation of the new threshold N cr mimics that of N cr , with the factor 
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K taken into account: 



e 

i=i 

~ M^r(|)c(|) 7 
Vcr (Kr(|)c(|)) 3 / 5 cr ' 

where N > N cr + A. 

Consider the following auxiliary problem: we wish to put some balls into the boxes 
Uj, j = 1, 2, ... , of multiplicities qj, leaving the box Uq aside. The overall number of balls 
is not specified in advance, and we should only observe the condition 

oo 

Y,JN 3 <M. (2.19) 

3=1 

Theorem 10 in [23] and Theorem 1 in [23] claim that in this problem the sum of all Nj is 

in most cases close to N cr . More precisely, one has the estimate 



p( N cr -f^Nj >A]< C m N~ 

V 7=1 ' 



N~ m (2.20) 



with some constants C m , m = 1, 2, 

Let G(L) be the number of ways to put exactly L balls into the boxes Uj, j = 1, 2, . . . , 
so that condition (12.191) is satisfied. Note that G(L) = for L > M, because 

oo oo 

/- (2.21) 

3=1 3=1 

Then the estimate (12.201) can be rewritten as 

Let J\f be the total number of ways to put N balls into the boxes Uq, U\, . . . with condition 
(12.41) being satisfied, and let jV(A) be the number of only those ways for which, in addition, 

\N -(N-N cr )\ > A. (2.23) 

One obviously has 

AT' 

M = G(L)F(N - L), (2.24) 

L=0 

where F(x) is the number of ways to put x balls into the box Uq of multiplicity q and 
N' = min{iV, M}. In a similar way, 

7V(A) = G(L)F(N — L). (2.25) 

0<L<N' 
\L-N cr \>A 

1 8 



Note that F(x) is a monotone increasing function. Hence we can estimate 
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Af > F(N - N') G ( L ) > 2 F ( N - N ') ( 2 ' 26 ) 

L=0 " L=0 

(The last inequality follows from (I2.22p and (ii).) Next, 

JV(A) < F(N) ^ F ( N ) ^ 2 - 27 ) 

0<L<N' |L-7V cr |>A 
|L-7V cr |>A 

By dividing fl2T27D by (E2BD , we obtain 

AT " F(N-N') £ L G(L) " m cr F(N-N') 1 j 
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in view of (12.221) . It remains to note that F(x) ~ Cx qo 1 with some constant C > 0, and 
hence 

F(N) / AT Nl/2 

F(N - AT' 

By substituting this into (12.281) . we obtain the desired estimate. The proof of the propo- 
sition is complete. 

At K > 1 the chemical potential /x in ( 12. 17ft is strictly less than zero, hence the 
compressibility factor will be greater than the value of (12.16!) . 
First, consider the graphs in Figs. [Hand E]for argon. 

If the vapor is saturated, then, at low temperatures, the number of clusters (dimers, 
trimers) rule, large. This decreases the total number of particles in the volume 

and increases the chemical potential, and hence the compressibility factor Z = PV sp /kT, 
P - is the pressure, V^ p is the specific volume, is increased. As the temperature increases, 
the number of clusters decreases and, at a certain temperature, the fraction of dimers 
becomes less than 7% (the Calo criterion). Then the compressibility can drop to 0.53. 

But since the saturated gas is in equilibrium with the liquid, the dimension can then 
decrease rather steeply and the compressibility factor (e.g., for argon) can decrease down 
to 0.25. It means that as the pressure increases, interaction takes effect. 

Thus, the formation of nanostructures in the other phase (the liquid one) plays a 
significant role, just as the formation of clusters in a gas. 

Let us now pass to the case of a constant temperature (Fig. [3j) [16] . 

We can assume that, instead of E, V 2 ^ 3 tends to infinity. And hence, in all the 
formulas of Bose-Einstein type from [2], [3], we can assume that b~ x ~ where T is 
the temperature, k is the Boltzmann constant, and d is the dimension. For the number b 
in [T71 [T8] to be dimensionless, let us introduce the effective radius a of the gas molecule 
(see below). 



3 Taking into account the pair interactions between 
particles 

Now consider the Hougen-Watson diagram given in Brushtein's textbook "Molecular 
Physics" [25]. The diagram reflects the dependence of the compressibility factor Z = 
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Figure 1: Thermodynamic properties of saturated argon. Z is the compressibility factor, 
Z = PV/kT; T is the temperature in Kelvin degrees. 
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Figure 2: Thermodynamic properties of saturated argon. P is the pressure in pascals, T 
is the temperature in Kelvin degrees. 
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Figure 3: T r = T/T c , and P r = P/P c are reduced temperature and pressure, respectively. 
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Figure 4: The original Hougen- Watson diagram. 

PV/ kTN on the pressure for different temperatures and was constructed by Hougen and 
Watson for seven gases: H 2 , N 2 , CO, NH 3 , CH 4 , C 3 H 8 , C 5 H 12 . Although the textbook 
states that attraction decreases compressibility, this, however, is obtained for Van-der- 
Waals gas under the condition |1 — Z\ -C 1, i.e., time as a compressibility factor decreases 
down to 0.2. 

From our understanding of ideal gas, it follows from the distribution (12. ip that Z can 
attain the value of 0.523 and, further, the compressibility factor must decrease only at 
the expense of the interaction. 

Phenomenological thermodynamics is based on the concept of pair interaction. More- 
over, it is implicitly assumed that 

there exists some one-particle distribution characterizing the field, to which all the 
particles contribute. 

They are interrelated. Formulas for the distribution corresponding to this mean field 
were rigorously obtained by the author in [26]. The equation that relates the potential of 
the mean field to pair interactions is called the equation of self-consistent (or mean) field. 
For the interaction potential $, it is of the form 

u(x,p) =u (x,p) + N J <$>{x - x\p-p') e{pV2m+u(x } p , ))(kT) ^ _ i dx'dp'. (3.1) 

This equation was rigorously justified only in the case of long-range interaction, in par- 
ticular, in [27] . 

In the case of a gas occupying the volume V and not subject to the action of external 
forces, Uq(x,p) = for x lying inside the volume V and Uq(x,p) = oo on the boundary of 
this volume. 

In what follows, we shall study only this case. For the interaction potential we take the 
Lennard- Jones potential or, for a greater coincidence with the experiment, the following 
potential: 

Hr) ^('l)""'- 6, (El.^) (3.2) 
v ; n-6 \6 J \r n r e J K J 

containing one more parameter n > 6. 

In the zeroth approximation, as o — > 0, the integral of the potential (13.11) with respect 
to r from some ro to oo. This integral substantially depends on ro- How must we choose 
r ? 
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Consider the scattering of one particle by another. Suppose that, as t — > — oo, the 
velocities of the two colliding particles are equal to v™ and v| n , respectively. This means 
that, at t — > — oo the trajectories of the particles approach straight lines. In terms of the 
variable r = r 2 — r% as t — > — oo the radius vector of the r-point asymptotically approaches 



the function r m = p + v m t, where pv m = and v 8 



Vj". The constant vector p is 



referred to as the target parameter. The quantity p is equal to the distance between the 
straight lines along which the particles would move if no interaction was present. After 
the collision as t — ► oo, the velocities of the particles are equal to v™' and vf* 1 *. This 
means that the radius vector r(i) asymptotically approaches the function r out = c + \ out t. 
The trajectories r m (t) and r out (t), which are straight lines, are said to be the incoming 
and outgoing asymptotes. The value of the relative speed in the in- and out-states is 
preserved; namely, |v m | = |v ou '| = v. 

The condition on the turning point r is of the form 
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Let is find the value of $(r) for the potential (13. 1| (this value depends on r ) 



$(0) 



$(|r|)dr = As I ( "— - "— ) dxdydz = / ( - ^ ) ^ 



ro 



poo 






/ 

'ro 




r 6 ) 



r/<r 



16tt 



-£<7 



9 V^o 



+ 




16 



3 Tq 9 7"n 



(3.4) 



where r/a = £ s . 

In particular, for r = a, the area is equal to — (8/9)ecr 3 a, where a = (4/3)p. For 
r = the area is equal to zero. 

In the nanotube, the target parameter p can be assumed to be zero. In this case, 
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(3.5) 



for p = 0, r = a. 

The Lennard- Jones potential can now be represented in the form $(r , r) = $(r (p), r), 
and hence, since \p\ = \p$ —Pj\, it follows that the equation for the dressed potential looks 
as follows: 



u(p, x) = N $(r (p - 77), \x - ^|)- 
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(3.6) 



Since the external potential is absent and the distribution depends on x in terms of the 
dressed potential only, we can assume that u(p,x) = u(p) does not depend on x. 



Making the change x — £ = y and integrating with respect to y from r to oo, we 
obtain 



u(p) 
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Here p <0 stands for the chemical potential. 

It should be noted that the probability of the event in which the particle x 2 occurs on 
the interval from r to infinity is not constant. It is obviously proportional to the time 
during which the particle is kept within the interval (x' 2 ,x 2 ), and this time is inversely 
proportional to the speed of x 2 with respect to the particle x±. One can readily see that 
this probability is equal to 
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Since the integral of the expression 
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(first with respect to x 2 and then with respect to dp) is equal to one, it follows that the 
probabilities are independent, and the distribution with respect to x 2 and p 2 is equal to 
the product of the distributions with respect to x 2 and to p 2 . 

Let us find the energy level below which the "condensate" appears. 

As is well known, the "turning point" r , the energy E = m(vi — t> 2 ) 2 , where V\ and 
v 2 are the velocities of two interacting particles, and the impact parameter p are related 
by 

Ep 2 



E 



$(r ) = 0. 



(3.9) 



1. The potential has the form — a/r A . Then the turning point r in the scattering 
problem is defined by the relation 



1-^ + 



a 
r^E 



0, 



where E = (pi — p 2 ) 2 /m is the energy of the particles an infinite distance apart and p is 
the impact parameter. Hence 
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~E ' 



and the solution is only possible if E is bounded below: E > 4a /p 4 . 
For E = 4a /b 4 , we have the expression 
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Figure 5: The Hougen- Watson diagram for nitrogen; c is the critical temperature. 



where r > p 2 /2. 

2. Suppose that the attraction potential is of the form 



$( r ) = -AU 
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EL 



min E = 4C/oc 6 max 



ro<p r. 
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Then we obtain 
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4 p 6 
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Note that, for E < -Emm, if there is no term corresponding to repulsion, then the "falling 
on the center" phenomenon occurs, i.e., the binding-together of the particles. 

3. For the function 
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$ = 4f/n — -^77 ^-1 



12 
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the equation for extremum points is of the form 

,3x 8 - 2p 2 x 6 - 6a 6 x 2 + 5a 6 p 2 



&(x) = 8ar e 
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The graph of $(r ) for the given impact parameter p = 2a is shown in Fig. [51 
For E less than some value of E , there appears a barrier whose the depth is -Emax — -^min 
and for which the probabilities of penetration of particles into the well PQ, due to 
thermal noise at a given temperature are well known. This implies that, below this 
barrier, the distributions that were obtained earlier for an ideal gas are false, because 
there exists a probability of penetration through the barrier at a given temperature. 
Thus, the maximum of the barrier is the natural minimum for the energy (pi —p2) 2 / m i 
below which we cannot use the distribution given above. 

In pp, [28], it was calculated how much time a particle stays in a well of height h and 
depth 5 at a given temperature (thermal noise). If the number of particles N tends to 
infinity, then this time is proportional to the number of particles occupying the well. It 
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is a strange fact that this is in agreement with our estimates [ 1 1 J of penetration into the 
condensate and yields the value E min (p) for the condensate; for a given impact parame- 
ter p, this value corresponds to the same turning point tq = 1.21 for the Lennard- Jones 
potential with n = 12, to r = 1.02 with n = 18, and to r = 1.28 with n = 7. 

But since p 3 for an ideal gas is of the order of 2V sp ^> p 3 , then, accordingly, we can also 
define E min from Eq. (13. 9p for prescribed values of r and p = ^/2V sp ; E min corresponds 
to the gases for which the interaction between particles is best described by the Lennard- 
Jones potential with a given n. If this is known, then E m[n can be determined in terms 
of V sp . 

Using V p ((pi — p2) 2 /m), we define the one-dimensional distribution of the difference 
of the momenta in the scattering problem [14J, [29] 



Vr, 



(Pi - P2f 



m 



1 



e {pi-p 2 ) 2 /mkT-n'/kT _ I 



e i 2 /mkT-ii> /kT 



(3.10) 



Here p min = ^mE min , and, as pointed out above, E min is defined for the scattering 
problem by an interaction in the form of the Lennard- Jones potential. 

The distribution (13.11) contains the chemical potential pi which is related to the chem- 
ical potential p for the distribution with a "dressed" potential by relation (13.111) below, 
which expresses the fact that, for a fixed scattering parameter p, the number of particles 
in the one-dimensional scattering problem is of the order of y/ N/2, where N is the total 
number of particles outside the condensate. 

The dressed potential u depends on the three-dimensional momentum p and is inde- 
pendent of the coordinates under the reduction to the scattering problem [22]. Therefore, 
the chemical potential pi is connected to the chemical potential p = p(p) of the problem 
on the distribution with a dressed potential by the relation 

3 



-7T 



p dp 



,(p2/2m+u(p,p)- l i)/kT _ I 



dp 



e p 2 /rnkT-fix/kT _ J 



(3-11) 



If u(p) is positive, then p min = 0, and hence pi < 0. 

In the integral equation of the mean field (I3.ip . we can drop the external potential Uq, 
because it is zero inside the volume, and finally obtain 



u(p,p) + C(p,p) = -vr 
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where V is the volume, r 
As presented in [26J, 



^pi/m-p^/kT _ l 

a 

r (p, E), and E = (p — p') 2 /m. 
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(3.12) 



p 2 dp 



e p 2 /2m+u(p,p)-p,/kT 

(3.13) 



When the scattering problem is considered in the whole space, then the distribution 
over the scattering section is uniform. But we restrict the problem by the volume V^ p . 
Then there is no uniformity due to the boundary, at least, outside the domain, where 

dZ Z 
dV sp V sp ' 

Under different assumptions, the distribution T(p) over p can be different (see j3U] (the 
Bertrand paradox), [3T] , [32] )• Moreover, Z is described by an expression of type (13.131) . 
averaged with respect to the distribution V(p) of the lines p apart in the ball of radius 



3kTV^J 3 JO J I 2m ^' K7 J e p 2 /2m+u( P ,p)-fi/kT _ ]_ eP 2 /2m+u( P ,p)-n/kT _ ]_ 



(3.14) 

where p = p(p). Therefore, Z can be taken at some mean point p me£lIl (T, V sp ). Then 

dZ dZ rfpmean 



dV sp dp mean dV sp 

and hence, 

dZ 

= — oo. 



dpmem 

If, at this point, the asymptotics as p — > of u(p, p mea n) is of the form —p 2 /2m + 
«(Pmean)| hip|, where da/dp mea _ n > 0, then this leads to the domain in which dZ/dp mea _ n = 
— oo, defining the A-transition to the condensate state and to the law 

dZ_ = Vs P 
dP kT 

for P > P\ and V sp > V\ (P\ and V\ depend on T). 

The equation for the dressed potential is of the form (we have omitted the chemical 
potential for simplicity) 

uM = r r«p-M-f)r^, , r r iv y _ (315) 



where 
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(3.16) 



0((p — ?]) 2 ) is nonzero only in the domain |p — p'\ 2 /m > E min , p 3 = V^ p , and C = C(p,p). 
Let us rewrite this equation in the form 

, , 1 f 00 F((p-r,)*)G{(p-r,)*)T?dT, {T dp y A 

U \PiP) 2 / e ( V 2 /2m+u( V , P ))/kT _ I ) / eP 2 /mkT_l\ U 
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For example, A. M. Chebotarev proposed the following distribution of the impact parameter: P(p < 



r) = (1 - r 2 ) 3 / 2 
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and make the replacement (p — rj) 2 = £ 2 . 
Then 



1// 00 g(g)de ^(e 2 )^ 
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We search for conditions under which the solution of this equation, as p —>■ oo, is of the 
form 

where dc/dp > 0. 

First, note that, by virtue of proofs and estimates similar to those given in the the- 
orems , we p ut the upper limit of the integral over 77 equal to infinity, because E in 
( 11. lip . p.!3p . ( 11.151) (which is different from E = (jp — p') 2 /m in the scattering prob- 
lem) is large, while the integrand is rapidly decaying, and the difference between the 

limit \JlmE and oo is less than the given estimates. 

But since we are concerned with the asymptotics of the solution u(p, p) as p — » oo, it 
follows that, as p' — > oo, the integral over p' must sufficiently rapidly converge. Therefore, 
this remark must be taken into account only for some exotic family of solutions. 

After the replacement indicated above, we express the term rf as rf = (£ —p) 2 . But 
since the function is symmetric, it follows that the integrati on of 2p£ over £ yields zero. 
Thus, we find that the term on the right-hand side of Eq. (I3.15P is proportional to p 2 . 
Now it suffices to equate the integral over £ as p — > oo to —p 2 /2m. Moreover, the choice 
of the constant c(p mean ) remains arbitrary. 

After the replacement, we obtain 

p 2 

u(p,p) = -— +w[p,p); 
2m 

here, as p — > oo, we have w(p, p) = c(p) | lnp| , and we can write an equation for the function 
w(p,p) > 0, dc(p meSun ) / dp mean > 0. Moreover, the phase A-transition, just as the minimal 
point of the condensate, depends on the power of the repulsive term in the Lennard- Jones 

potential as well on the quotients 7 = a/Vjp^ 3 and a = U /kT. The points 7 crit and the 
minimal point a, corresponding to dZ/dV sp = —00 are called X-critical. As the pressure 
w(p,p) increases above the point of the A-transition, V^ p remains unchanged. This implies 
that the volume V decreases as the pressure increases, but, simultaneously, the number 
of particles outside the Bose condensate also decreases. It is possible that all the particles 
became dimers. and hence the total number of particles has decreased. Further, they all 
became trimers, etc. The volume V has decreased, while the specific volume V sp remained 
constant — this is the law of the Bose condensate for classical gases or, more precisely, is 
the law of cluster formation. 

The energy of the A-point of the logarithmical form appears as T = T cr . As one can 
see in Fig. [5] and by virtue of 



dp . d 2 p 



\T=Tci,p=pci A\ri \T=T CI ,p=p C t 







the coefficient of incompressibility k = — (din V sp )/dp turns to infinity, and the compress- 
ibility factor decreases steeply. This arouses a wave of compressibility (a shock wave), and 
therefore an additional term c\p\, where c is the speed of sound, must o ccur in energy. 

This term for u(p, p) slows down the decrease of Z, and equation (13.251) eliminates 
the shock wave as well as this term. At that instant for the derivative of heat capacity 
the so-called A-point occurs. This effect is shown more obviously in Fig. [31 
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This phenomenon, as well as consequences of the Pontryagin-Andronov-Vitt theorem, 
does not follow from classical mechanics but occurs when noise and fluctuations are taken 
into account. Therefore it does not follow from formulas for the dressed potential, although 
the equations of "collective oscillations," as well as "equations of variations," are related 
to the dressed potential (33J, [M] . 

To illustrate this we use both wave and quantum equations. The wave equation of 
sound propagation has the form 

whereas the Schrodinger equation is 

f h 2 1 

ih— = < -— A + u{x) (3.19) 



or in the iterated form 



dt 1 2m 



2 d 2 ^ f h 2 



2 



As follows from formula (25) [31] and [35] , 

/ l 2 ^ = /iVA$-f-^A + H (x)l ^ + 6(/i 2 )^, (3.21) 
ot z l 2m J 

where 0(h 2 ) is an operator such that Oe s(yX,t ^ h = 0{h 2 ) for any C°°-smooth S(x,t). 
For u(x) = and ^ = e ( l / h )(p x - Et ) this implies 

E 2 = c 2 p 2 + ^- 2 , (3.22) 

which coincides with the spectrum obtained by N.N. Bogolyubov for the weakly nonideal 
classical gas [37] . This author has established in 1995 that this spectrum has quasiclassical 
rather than quantum nature [33j [3l] and, based on a dependence from the capillary 
radius derived in [2H], applied these results to the classical gas in nanotubes [3H] • These 
predictions are justified by autho ritative experimental data [^0] . 

The phase-space boxes HI . 23[) are chosen to be invariant with respect to the Hamilto- 
nian system corresponding to the Hamiltonian function H (p, q) . In the context of quantum 
chaos |2~Tj I^Tj, t his corresponds to a kind of generalized ergodicity. Moreover, it follows 
from IMp 55 ] ES] that the Hamiltonian 

V<V + H 2 (p,q) (3.23) 
corresponds to the Arnol'd diffusion in a self-consistent field. 



In this case, the equation for the dressed potential has the form 



u(p,p) + C(p,p) = -vr — 
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exp(((p — p') 2 /m — pi)/kT) — 1 
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dp 



Pmi „ exp((p 2 /m)/A;T) - 1 
p 2 dp 



exp 

((^cp' 2 + p' 4 /4m 2 + u(p, p))/fcT) - 1 



(3.24) 



, (3.25) 



where p and p are given, C (p, p) is a constant depending on p and p, V^ p is the specific 
volume, r = r (p, J5), and p min is determined by condition (13. lip for p x = p — 0. 

After the change of variable w(p, p) = u(p, p mea n) — c|p| we find the values of T and 
V sp for which w(p,p mcan ) w 0(p 3 ) as p(p mcan ) -> 0. 

It follows from the above theorems that as iV — > oo, it is necessary to introduce a 
parameter k in the exponential in (13.1 ip in the left-hand side and a parameter K\ in the 
exponential in the right-hand side. Since p m { n — > 0, we have K\ 3> K, and hence the kernel 
of the integral operator tends to the 5-functoin as k — >• 0: 8{p — p' + p m m)- To cancel the 

term — in the right-hand side as p — > 0, it is necessary to satisfy the relation between 

Pmm and p m ean- 

Then the leading term of the T derivative of TZ, which contains the heat capacity C v , 
will feature a logarithmic dependence characteristic for a A point. Indeed, as p(p me an) — ► 
we have 

Q(p 3 ,Pmean)p 2 ^P f dp . . 

( e O(p3, Pmcan ) _ ^ 2 ~ J p ■ ^-^f 

Observe that although the quantum equations for the self-consistent field go over into 
the classical ones as h — > 0, the equations of variations for the quantum mechanical 
equations of the self-consistent field assume in the same limit an extra term with respect 
to the classical equations of variations ( the eq uations of collective oscillations, see [39l 
1.1]). This gives rise to the Hamiltonian (13.231) . 

One can assume that at temperatures below the A point ergodicity turns over into a 
KAM situation, making supefluidity of a classical gas possible in a very thin nanotube 
capillary. Thus the temperature of the A point can be regarded as the crossover point 
between the generalized ergodicity and the KAM dynamics. 
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